Measurements of scattering observables for the pd break-up reaction 
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Abstract. High-precision measurements of the scattering observables such as cross sections and analyzing 
powers for the proton-deuteron elastic and break-up reactions have been performed at KVI in the last two decades 
and elsewhere to investigate various aspects of the three-nucleon force (3NF) effects simultaneously. In 2006 an 
experiment was performed to study these effects in p + d break-up reaction at 1 35 MeV with the detection system, 
Big Instrument for Nuclear polarization Analysis, BINA. BINA covers almost the entire kinematical phase space 
of the break-up reaction. The results are interpreted with the help of state-of-the-art Faddeev calculations and are 
partly presented in this contribution. 



In the last few decades, several nucleon-nucleon poten- 
tials (NNPs) have been studied extensively to describe the 
properties of bound nuclear systems by comparing high- 
precision two-nucleon scattering data with modern poten- 
tials based on the exchange of bosons II 11213 141 . The mod- 
ern NNPs reproduce the world database with a reduced^ 2 
close to one and have, therefore, been accepted as high- 
quality potentials. Deficiencies of theoretical predictions 
based on pair-wise nucleon-nucleon potentials have been 
observed in three-nucleon scattering observables. For ex- 
ample, exact solutions of the Lippmann-Schwinger equa- 
tions (Faddeev calculations) |5| solely based on modern 
NN interactions fail to describe high-precision differential 
cross sections of proton-deuteron elastic scattering at in- 
termediate energies obtained at many laboratories includ- 
ing KVI [6 7 8 9 1, Research Center for Nuclear Physics 
(RIKEN) |[T4lT5l and RCNP Q6). Calculations based on 
NNPs including 27r-exchange type three-nucleon forces (3NFs) 
remove this discrepancy for a large part and lead to a good 
description of the measured cross sections for energies be- 
low 100 MeV/nucleon. However, the description of spin 
observables such as vector and in particular the tensor an- 
alyzing powers is not satisfactory and an inclusion of 3NFs 
is not sufficient to remedy the observed discrepancies 111 1I6I7I8I . 

The break-up reaction has a rich phase space which al- 
lows a systematic study of the 3NF. Predictions show that 
large 3NF effects can be expected at specific kinematical 
regions in the break-up reaction. For other phase space re- 
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gions, the effect can be small, which make them, therefore, 
suitable for benchmark studies. In elastic and break-up re- 
actions, precision data for a large energy interval for the 
differential cross section and analyzing power have come 
from recent experimental studies at KVI 1617181 10191 121 131 . 
Break-up cross sections and analyzing powers for a beam 
energy of E d lab = 130 MeV have been published in Refs. II 1 71 1 81 1 9120 
For a further study of 3NF effects, we have measured the 
break-up cross sections and vector analyzing powers for a 
proton beam energy of E p l b = 135 MeV. In this experiment, 
cross sections and vector analyzing powers in p + d — > 
p + p + n reaction were measured using a polarized beam 
on a liquid-deuterium target ll22l . The results are compared 
with predictions derived from state-of-the-art Faddeev cal- 
culations and are partly reported here. 

The proton-deuteron break-up experiment described in 
this paper was performed using the recently-developed de- 
tection system, BINA. The energy correlation between the 
two outgoing protons, E\ vs. Ez, as shown in Fig. Q] was 
studied for a large number of kinematic configurations. To 
obtain the cross section and the analyzing power, several 
slices along the S -curve were made. The angular ranges 
for the integration of events was chosen to be AO\-A02-^° 
and A(p\2-\Q°, which were wide enough to have sufficient 
statistics, while variations of the cross section within these 
ranges are small. In this way, the experimental cross sec- 
tions can be directly compared with the theoretical predic- 
tions calculated for the central values of a specific config- 
uration. The projection of the indicated region in Fig.Q]on 
the line perpendicular to the S -curve is shown in Fig. [2] 
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tribution is fitted using a third-order polynomial together 
with the main peak represented by a Gaussian function. 
The number of events underneath the Gaussian peak will 
be referred to as the number of break-up events and will 
be used to obtain the cross section and analyzing power 
after making efficiency corrections such as the one from 
hadronic interaction. 

Figure [3] shows the cross section, dQ f d ^ ds [//b/(sr 2 MeV)], 
as a function of S [MeV] for the symmetric and coplanar 
configuration with (0 u O 2 ,(f>n) = (25°, 25°, 180°). In this 
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Fig. 1 : The coincidence spectrum between the ener- 
gies E\ versus E2 of the two protons registered at 
(0 1 ,6/ 2 ,012)=(25 O ±2 O ,25 O ±2 O ,18O O ±5°). The solid line shows 
the kinematical curve, the so-called 5 -curve calculated for the 
central values of the experimental angular ranges. 



This spectrum contains break-up events around channel 
and background from accidental events and contributions 
resulting from the hadronic interaction of particles in the 
scintillator material of the detector. This projection axis is 
denoted the D-axis. The crossing point of the S -curve with 
the D-axis defines the zero point. 
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Fig. 2: The projection of the slice chosen in Fig. [T] along the 
D-axis. The solid line is the sum of a third-order polynomial 
background function and a Gaussian distribution representing the 
peak. 



The accidental background is located at channels higher 
than the main peak in Fig. [2] and is already subtracted ex- 
ploiting the relative time-of-flight information between the 
two protons. The contributions from hadronic interaction 
are located at lower channels on the D-axis which corre- 
spond to events below the S -curve in Fig. Q] These events 
are mostly "true" break-up events. At this stage, this con- 
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Fig. 3: The cross section is plotted as a function of S [MeV] 
for the kinematical configuration, (61, 82, <f>n) - (25°, 25°, 180°). 
Lines represent Faddeev calculations from the Hannover-Lisbon, 
and Bochum-Krakow groups. The dotted line represents the cross 
section using the CD-Bonn two-nucleon potential, the solid line 
shows the CD-Bonn+TM' calculation. The dashed line repre- 
sents the results of a calculation by the Hannover-Lisbon group, 
which is based on the extended CD-Bonn potential, including a 
virtual A excitation in a coupled-channel approach. 



figure, the lines represent the predictions by the Bochum- 
Krakow and Hannover-Lisbon theory groups [23 24 25 26 .27]. 
The dotted line is the result of calculations using the CD- 
Bonn two-nucleon potential and the solid line presents the 
calculation including the three-body force, TM', as well. 
The dashed line represents the result of a calculation by the 
Hannover-Lisbon group, which is based on the extended 
CD-Bonn potential including a virtual A excitation in a 
coupled-channel approach. 

In the following, we describe the determination of the 
vector analyzing power, A y , by using a polarized proton 
beam and by measuring the induced asymmetry in the cross 
section. The relation between dcr s - dcr 1 ^ and the unpo- 
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larzied cross section, d<x°, is 



da tl = dcr°(l + p tl ■ Ay ■ cos <f>), 



(1) 



for an incoming proton beam with spin up (T) or spin down 
(I) and a vector analyzing power, A y . Here <p is the angle 
between quantization axis for the beam polarization and 
the normal to the scattering plane in the laboratory frame 
of reference. From the two cross sections, dcr T and dcr 1 , 
with polarizations, p T and p\ the analyzing power can be 
obtained from the asymmetry, 



Ay COS (j) = 



da^ - dcr^- 
dcrtpl - dcr^p 



(2) 



The reaction asymmetry can be measured by exploiting 
the distribution of events obtained with beam polarizations 
up and down together with the values of the beam polar- 
ization in these two modes. This gives a periodic function 
in cp with an amplitude that corresponds to A y . Figure [4] 
shows the asymmetry as a function of (f> for a particular bin 
in S . By exploiting the asymmetry distribution for each 
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Fig. 4: The asymmetry, AyCoscf), in the break-up reaction as a 
function of the azimuthal scattering angle of one of the protons, 



S -bin, the vector-analyzing power, A y , is obtained for ev- 
ery kinematical configuration, (61,62, <f>\2)- Figure ^repre- 
sents the analyzing -power for the configuration {61,62, <p\i) 
= (25°,25°,180°). The various lines represent calculations 
from the Hannover-Lisbon and Bochum-Krakow theory gr- 
oups as explained before. 

We determined the cross sections and analyzing pow- 
ers for configurations in which 14° < #12 < 30°, and the 
azimuthal opening angle, <pi2, varied from 20° to 180° in 
steps of 20°. Fig. [6] shows the results of the fixed combina- 
tion (81,62) = (25°, 25°) with different azimuthal opening 
angles, <p\2, as a function of S . The top panels depict the 
cross sections and the bottom panels show the analyzing 
powers. 

The presented error bars in all the figures are statistical. 
The main source of systematic uncertainties for the break- 
up cross sections are: the uncertainty in the target thickness 



Fig. 5: Analyzing powers of the break up reaction for the coplanar 
kinematics, (6\, d 2 , <p\%) = (25°, 25°, 180°), as a function of S . For 
the explanation of the curves, see Fig. [3] 



(3.85+0.2 mm— »~5%), the correction for the hadronic re- 
action efficiency for both protons obtained via a GEANT-3 
simulation (92+3%— > ~6% for two protons), the correc- 
tion for the efficiency of the MWPC (92+1 % -> 2% for two 
protons), and the correction for the geometrical inefficien- 
cies obtained via GEANT-3 simulations which is at most 
12% for small azimuthal opening angles, tpu = 20°, and 
at most 2% for the larger azimuthal opening angles. Alto- 
gether, by adding the systematic uncertainties in quadra- 
ture, the maximum systematic uncertainty for cross sec- 
tions at small azimuthal opening angles is less than 14% 
and at larger azimuthal opening angles less than 9%. The 
systematic error for the analyzing power stems primarily 
from the uncertainty in the measurement of the beam polar- 
ization via the proton-deuteron elastic-scattering reaction. 
For instance, the beam polarization for the down-mode has 
been measured at a value of ~0.70 ± 0.04, which gives rise 
to a 6% systematic uncertainty in the analyzing power. 

The predictions of the Faddeev calculations using dif- 
ferent NN and 3NF models are added to every panel with 
different line colors and styles. The blue (long-dashed), red 
(dash-dotted), green (dashed), and black (solid) lines cor- 
respond to calculations based on CDB (NN), CD B+TM' 
(3NF) from the Bochum-Krakow group I23I24I25I . CDB+A 
(3NF), and CDB+zl+Coulomb calculations from the Hann- 
over-Lisbon group [26 27], respectively. Here, all theory 
curves have been calculated in a fully non-relativistic frame- 
work with non-relativistic observables and, therefore, the 
length of S for these calculations is slightly shorter than 
that in the data. The typical difference in length of S for 
the relativistic and non-relativistic kinematics is less than 
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Fig. 6: The cross sections and the analyzing powers at (9 l ,8 2 ) = (25°, 25°) as a function of S for different azimuthal opening angles. The 
error bars reflect only statistical uncertainties. The blue (long-dashed), red (dash-dotted), green (dashed), and black (solid) lines show 
predictions of Faddeev calculations using CDB (NN), CDB+TM' (3NF), from the Bochum-Krakow group I23I24I25I . CDB+zf (3NF) 
and CDB+zf+Coulomb calculations from the Hannover-Lisbon group |26 27|, respectively. 
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1-2 MeV, depending on the azimuthal opening angle, <pi2- 
This difference is less than the experimental resolution in 
S of 4 MeV (FWHM), and we, therefore, did not transform 
the non-relativistic S -curves to the relativistic ones. 

For the configurations at large azimuthal opening an- 
gle, > 40°, and taking the systematic uncertainties into 
account, a reasonable agreement is observed between the 
cross-section data and the corresponding theoretical pre- 
dictions. For the configuration with a small relative az- 
imuthal angle, the picture changes. Here, the measured cro- 
ss sections show a large discrepancy with a calculation 
which includes the TM' 3NP. In this region, the CDB+zl+C- 
oulomb calculation has a smaller deficiency when com- 
pared with the experimental data, but the deficiency is still 
large for small value of polar angle, <p\2 = 20°, as shown 
in Fig. [6] 

For the analyzing powers, the major discrepancies be- 
tween the data and the theoretical calculations arise at small 
azimuthal opening angles. In this range, the predictions 
based solely on a NN potential are closest to the data, al- 
though, the disagreement is still significant. The inclusion 
of 3NPs increases the gap between data and predictions as 
can be seen from bottom panels of Fig. [6] The contribu- 
tion of the TM' 3NP appears to be larger than the implicit 
inclusion of the A resonance by the Hannover-Lisbon the- 
ory group. It is interesting to note that a similar, but even 
larger, discrepancy has been observed in a break-up study 
at an incident beam energy of 190 MeV [ 10 1. 

This paper discusses some of the preliminary results 
of a proton-deuteron break-up experiment carried out with 
an incident proton beam of 135 MeV. The data were taken 
using a nearly An detection system, BINA, and exploiting 
a polarized beam of protons. With this, precision differ- 
ential cross sections and analyzing powers were measured 
and compared to Faddeev calculations based on modern 
two- and three-nucleon potentials. The large coverage of 
the detection system provides an ideal tool to systemati- 
cally explore the rich phase space of the break-up reaction. 
We have identified various configurations at which signifi- 
cant discrepancies are observed between our data and pre- 
dictions by Faddeev calculations based upon state-of-the- 
art potentials. Intriguing deficiencies are observed for the 
analyzing power for configurations at which the relative 
energy between the two outgoing protons becomes small. 
The discrepancies cannot be explained by the Coulomb in- 
teraction and higher-order relativity, since these effects are 
accounted for in the present state-of-the-art calculations. 
Therefore, the data provide an ideal basis to develop a bet- 
ter understanding of three-nucleon force effects in few- 
nucleon interactions. 



7. K. Ermisch et al, Phys. Rev. C 68 051001(R) (2003). 

8. K. Ermisch et al, Phys. Rev. C 71 064004 (2005). 

9. H. Amir-Ahmadi et al., Phys. Rev. C 75 041001(R) 
(2007). 

10. H. Mardanpour, Ph.D. thesis, KVI, University of 
Groningen, 2008. 

11. R. Bieber et al, Phys. Rev. Lett. 84 606 (2000). 

12. A. Ramazani-Moghaddam-Arani et al., Phys. Rev. C 
78 014006 (2008). 

13. M. Eslami-Kalantari, Ph.D. thesis, KVI, University of 
Groningen, 2009. 

14. H. Sakai et al, Phys. Rev. Lett. 84 5288 (2000). 

15. K. Sekigushi et al, Phys. Rev. C 65 034003 (2002). 

16. K. Sekigushi et al, Phys. Rev. Lett. 95 162301 (2005). 

17. A. Biegun et al, Acta. Phys. Pol. B 37 213 (2003). 

18. St. Kistryn et al, Phys. Rev. C 68 054004 (2003). 

19. St. Kistryn et al, Phys. Rev. C 72 044006 (2005). 

20. St. Kistryn et al, Phys. Lett. B 23 641 (2006). 

21. E. Stephan et al, Phys. Rev. C 76 057001 (2005). 

22. N. Kalantar-Nayestanaki et al, Nucl. Instr. and Meth. 
in Phys. Res. A 417 215 (1998). 

23. H. Witala et al, Few-Body Systems 15 67 (1993). 

24. H. Witala et al, Phys. Rev. Lett. 81 1 183 (1998). 

25. H. Witala et al, Phys. Rev. C 63 024007 (2001). 

26. A. Deltuva et al, Phys. Rev. C 71 054005 (2005). 

27. A. Deltuva et al, Phys. Rev. C 72 054004 (2005). 



References 

1. V.G.J. Stoks et al, Phys. Rev. C 48 792 (1993). 

2. V.G.J. Stoks et al, Phys. Rev. C 49 2950 (1994). 

3. R.B. Wiringa et al, Phys. Rev. C 51 38 (1995). 

4. R. Machleidt, Phys. Rev. C 63 024001 (2001). 

5. W. Glockle et al, Phys. Rep. 274 107 (1996). 

6. K. Ermisch et al, Phys. Rev. Lett. 86 5862 (2001). 



